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Abstract 

Finite-temperature spectra of heavy quarkonia are calculated by combining potential 
model and thermofield dynamics formalisms. The mass spectra of the heavy quarkonia 
with various quark contents are calculated. It is found that binding mass of the quarko- 
nium decreases as temperature increases. a 

Quarkonia; finite-temperature spectra; thernofield dynamics. 



1 Introduction 

Recent progress on creation hot and dense matter using LHC and RHIC lead to new challenges 
in particle physics. In particular, it has been reported recently that creation of quark-gluon 
plasma has been achieved in the relativistic heavy ion collisions experiments [1]. A key signal 
in the experiments on quark-gluon plasma formation is the J/\E' suppression [1, 2, 3, 4] which 
implies that the quarkonium properties are changed under the influence of in-medium effects, 
such as high density and temperature. Therefore precise knowledge or qualitative estimates of 
these changes is of importance for such experiments. 

Heavy quarkonia have been the subject of extensive theoretical and experimental research 
for last three decades [5]- [10]. It is well-known that the properties of mesons and nucleons 
change in media of high density and at finite temperature. Even the coupling constants of 



1 



meson-nucleon interaction are altered as the temperature increases to the deconfining tempera- 
ture. The present study is focusscd on the properties of quarkonia at finite temperature [11, 12]. 
Several models have been used to describe the spectra and other properties at zero temper- 
ature [9]. Among others, nonrelativistic QCD [13, 14, 15], the lattice QCD [14, 16, 17] and 
the potential nonrelativistic QCD [18] lead to considerable progress in understanding the na- 
ture of quark-quark and quark-antiquark interactions and predicting various properties such 
as decay, mass spectra, fragmentation etc. Also, a quantum mechanical approach based on 
potential models [5, 7, 8] can be used to obtain quarkonium spectra for various quark contents. 
In this case one may use Schrodinger [5, 7] or Bethe-Salpeter [9, 19] equations for calculation 
of quarkonium spectra. Despite the considerable progress made in quarkonium physics most 
treatments have been limited to considering zero-temperature case. There are few papers in 
the literature with finite-temperature treatment of quarkonium properties. 

However, in medium effects on quarkonium properties should be taken into account in con- 
sidering quarkonia in quark-gluon, hadronic and nuclear matter. In experiments on quark-gluon 
plasma formation through high energy heavy ion collisions, heavy quarkonium states are the 
simplest probes that allow us to test the structure of the QCD vacuum. Especially, here is 
of interest understanding of the behavior of heavy quark boTind states in a strongly interact- 
ing media close to deconfinement temperature. Also, strict treatment of the J/^-suppression 
due to color-screening should be also carried out with account of finite temperature. Few 
attempts have been made to use finite-temperature lattice QCD [20], perturbative QCD at 
finite temperature [21] and to use color screening models [22] in order to calculate quarkoni- 
um properties at finite temperature. Also, a hadron string model has been applied to find 
temperature-dependence of charmonium [2]. In addition, the prediction of the quarkonium 
spectra with account of finite-temperature effects is of importance for ongoing and forthcoming 
experimental studies on finding of hadrons in dense and hot matter. 

In this paper we present an approach for obtaining quarkonium spectra at finite temper- 
ature. Combining potential model approach and a real-time finite-temperature field theory, 
thermofield dynamics, we calculate the mass spectra of charmonium at finite temperature and 
compare it with corresponding zero-temperature ones. Finite-temperature mass spectra of the 
heavy quarkonia have been recently calculated on the basis of solution of the temperature- 
dependent Bethe-Salpeter equation [3] and using hadron string model [2]. 

This paper is organized as follows: in the next section we give basic descritpion of the Ther- 
mofield Dynamics. In section 3 thermofield dynamics formalism is applied for the calculation of 
the quarkonium energy spectra at finite temperature. Section 4 describes obtained numerical 
results and their discussion. 

2 Thermofield dynamics 

Thermofield dynamics (TFD) is a real time operator formalism of quantum field theory at finite 

temperature in which any physical system can be constructed from a temperature dependent 
vacuum |0(/3)) which is a pure state [23, 24, 25, 26]. The thermal average of any operator 
is equal to the expectation value between the pure vacuum state |0(/3)) defined by applying 
Bogolyubov transformations to the usual vacuum state. TFD is based on two features. The 
first one is the doubling of the Fock space such that the original Fock space and its double are 
defined non-tilde and tilde space respectively. All operators are also doubled and the finite- 
temperature creation and annihilation operators are constructed by Bogolyubov transformation 
between tilde and non-tilde operators. This is the same procedure in writing down the vacuum 



2 



state at finite temperature. Mathematically, the field operators have the following properties: 

{AiAj) = Aj^Aj, 
{cA-i -\- Aj) — c Ai -\- Ajj 

i^) = 
(Ai) = Ai, 

[M,] = 0. 

It is well known [27] that in an algebraic approach the doubled set of operators may be con- 
sidered as a set of operators that relate to the physical observables, O, and a second set that 
are generators of symmetries, O. The hat operators are responsible, in particular, for time 
development and are needed for the purpose of scattering, decay and any transitions between 
states. The physical observables lead to the quantities that are measured in experiment. Both 
for physical observables and generators of symmetry after Bogolyubov transformations leading 
to finite-temperature creation and annihilation operators and to a pure vacuum state only the 
non-tilde operators are required for getting the appropriate matrix elements. 

For a given Hamiltonian H{a, a^) written in terms of annihilation and creation operators, 
TFD prescription implies doubling, thus constructing a new Hamiltonian H{a, a^, d, a^) 

H{a, a\ a, d^) = H{a, a^) — H{d, d^). (1) 

Then applying Bogolyubov transformations, which are given by 

a — a{P)cosh9 + d^{P)sinh9 

— a\l3)cosh9 + d{(3)sinh6 
d = a\j3)sinh9 + d{j3)cosh9 

— a{(3)sinh9 + d\(3)cosh9 

where 

= sinh'9 ^ (e^ - iy\ 

the Hamiltonian, Eq. (1) can be written in a temperature-dependent form. Physical observ- 
ables to be obtained using this presciption are described in terms of non-tilde annihilation and 
creation operators. In the next section we apply TFD prescription for the calculation of the 
energy spectra of heavy quarkonia. 

3 Finite-temperature spectrum 

Within the framework of the potential model, heavy quarkonium is described by the Schrodinger 
equation which is given by 

where H is the quarkonium Hamiltonian: 

2 r 
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with Z — and ctg being strong coupling constant. 

Multiplying both sides of this equation by and introducing the following time scaling 



we have 



r 2^ X 3 2lT 

I— [-r^ agV + Xr-^ + Vor^]"^. 

OT 2 



Then using the substitution 

this equation can be written in time-independent form as: 

Separating angular and radial variables we get 

^'h^y^^ + \{l + 1) - + + For']* = 



(2) 



Introducing annihilation and creation operators as 

1 d 1 



a = —;=— H ^r, a 



t _ _ 1 ^ 1 



Eq. (2) can be written in terms of these operators: 

[-i(a-at)(a+at)2(a-at)--^(a+at) + ^(a+at)3 + i/(Z+l) + l^ (« + a^)"\^ = (3) 

Thus, the quarkonium Hamiltonian in terms of operators a and can be written as 

F(a, at) = -l(a-at)(a+at)2(a-at)-^(a+at) + ^(a+at)=^ + ^/(/+l) + ^^ (a + a^' • (4) 
8 V2 2v2 2 L ^ ' 

Constructing the following matrix on the harmonic oscillator basis as 

E(n,n') = (n' | F(a,at) | n) = 
= -| .Jn{n-l){n-2){n-3)5n,n'-4 + ^{n +l){n + 2){n + 3)(n + 4)5„,„.+4 - 

-(2n2 + 2n + 3 + 4/(/ + 



+ 



V8 



+3 (nVn5n,n'-i + {n + l)Vn+ 15„,n'+l) - ^ (\/ri5n,n'-i + \/riTT5n,n'+i) + 



+ 



Vb 



^n(n - l)5„,n'-2 + y'(n + l)(n + 2)5„,„/+2 + (2n + 1)6„ 



To find finite-temperature spectrum of the heavy quarkonium we should apply the TFD pre- 
scription to the Hamiltonian, Eq. (4). Then we have 



H = H-H, 



(6) 



and applying Bogolyubov transformations we get 

H = -l{AiA2A2Aicosh^e + BiB2B2Bisinh*e + [A1A2B2A1 + A1B2A2A1 - A1A2A2B1- 

o 

-BiA2A2Ai]cosh^esinhe + [B1A2B2B1 + B1B2A2B1 - B1B2B2A1 - AiB2B2Bi]coshesinh^e+ 
+[AiB2B2Ai - A1A2B2B1 - A1B2A2B1 + B1A2A2B1 - B1A2B2A1 - BiB2A2Ai]cosh^esinh^e}- 

^ [A2Cosh9 + B2sinh9] H p[^2AAcos/i^6' + B2B2B2sinh^0+ 

V 2 2v 2 

+(^2^252 + ^252^2 + B2A2A2)cosh^esinhe + {B2B2A2 + B2A2B2 + A2B2B2)coshesinh^e]- 
~Vq (^A2A2Cosh'^e + B2B2sinh'^9 + {A2B2 + B2A2)cosh9sinh9^ 

where 

= a(/5) - at (^p) ^ A2 = a{p) + (p) , 
Bi^a{/3)-d^{/3), B2 = h{f3) + h\/3). 
Again, the (finite-temperaure) energy spectrum is obtained by diagonahzation of the matrix 

E = (n', n' I H(e, a((3), a\(3), a((3), a\(3)) \n,n)^ 

where K,L,M,N,0 are defined as: 

K = 4l{l + l)Sn,n' + (cosh^e + sinh'^e) [{2n^ + 2n + 3)5„,„/- 

-yjn{n - l){n - 2){n - 3)5„,„/_4 + ^(n + 1) (n + 2) (n + 3) (n + 4)5„,„.+4 j , 

L = {cosh^9sinh9 + cosh9sinh^6^ ^4n^(5„_„/_i + 4(n + l)^(^n,n'+i^ ; 

M = 2cosh'^esinh'^e^{4n'^ + An - l)5n,n'+ 

+ (n(n - 1) - ^n(n-l)) 5„,„'_2 + (^{n + l)(n + 2) - ^J {n + l){n + 2)5n,n'+2^ , 
N = {cosh9 + sinh9){y/n6n^n'-i + + 15n,n'+i), 
O — {cosh^O + sinh^O) i^n^/n5n,n'-i + 3(n + l)\/n-\- l5„,„'+i+ 

l)(n-2)5„,„'_3 + ^(n + l)(n + 2)(n + 3)5„,„/+3^ , 

P = - l)5n,n'-2 + yj{n + l){n + 2)5n,n'+^ (cosh^O + sin/i^^) + 

+2 + (n + coshOsinhO. 

We can obtain finite-temperature energy eigenvalues by numerical diagonahzation of this 
matrix. 



4 Results and discussion 



In Table 1 we present the charmonium mass spectrum, calculated by numerical diagonalization 
of the matrix Eq. (5) for various temperatures. As is seen from this table, the masses of 
the charmonium states decrease by increasing temperature. This result is in correspondence 
with the result of recent calculations of finite-temperature mass spectrum of charmonium [2], 
where a smooth decreasing of the charmonium mass with increasing of temperature up to the 
deconfinement temperature is observed. 

Table 1: Mass spectra of charmonium (GeV) at finite temperature (GeV) in rric = 1.46 GeV, 
as = 0.39 {Z = fa,), A = 0.2 GeV^, Vo = -0.2 GeV. 





n 


T = 


T = 0.1 


T = 0.15 


T = 0. 


1 


(1=0) 


2.481 


2.447 


2.025 


1.966 


2 


(1=0) 


2.901 


2.839 


2.441 


2.353 


3 


(1=0) 


3.812 


3.515 


3.401 


3.175 


4 


(1=0) 


4.460 


4.290 


4.255 


4.167 


5 


(1=0) 


5.277 


5.264 


5.147 


5.098 


1 


(1=1) 


2.934 


2.888 


2.682 


2.566 


2 


(1=1) 


3.382 


3.290 


3.027 


2.964 


3 


(1=1) 


4.108 


3.908 


3.809 


3.665 


4 


(1=1) 


4.586 


4.484 


4.451 


4.255 


5 


(1=1) 


5.915 


5.562 


5.330 


5.315 


1 


(1=2) 


3.799 


3.750 


3.487 


3.350 


2 


(1=2) 


4.323 


4.219 


4.135 


3.986 


3 


(1=2) 


4.571 


4.523 


4.428 


4.256 


4 


(1=2) 


5.342 


5.272 


5.023 


4.919 


5 


(1=2) 


5.999 


5.884 


5.753 


5.628 



In Fig. 1 is the pictorial presentation of the mass spectrum for some of the charmonium 
states. As it is clearly seen from this figure for all presented states, the mass at T = 0.15 GeV 
is considerably smaller than it is for T = 0.1 GeV. However, the difference of masses between 
T = 0.2 GeV and T = 0.15 GeV is smaller. It should be noted that in our calculations we 
didn't take into account spin and color-spin effects. 

Thus we have treated charmonium spectra at finite temperature using the thermofield dy- 
namics formalism. Motivation for the treatment of quarkonium properties with account in 
medium effects such as finite-temperature and density comes from recent advances in quark- 
gluon matter creation experiments, where some states of charmonium are considered as a key 
signature of such matter. In addition, knowledge of the hadron spectra at finite-temperature 
is important for the study of hadrons in nuclear and hadronic matter. So, the obtained results 
can be useful in forthcoming LHC experiments on the quark-gluon plasma creation. 
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Figure 1: The mass spectrum of cc charmonium (GeV) at finite temperature (GeV). 
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